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Abstract: Hydatid cyst caused by Echinococcus granulosus is one of the most important parasitic diseases around the 
world and many countries in Asia, including Iran, are involved with this infection. This disease can cause high mortality in 
humans as well as economic losses in livestock. To date, several molecular methods have been used to determine the ge- 
netic diversity of £ granulosus. So far, identification of £ granulosus using real-time PCR fluorescence-based quantitative 
assays has not been studied worldwide, also in Iran. Therefore, the aim of this study was to investigate the genetic diversity 
of £ granulosus from center of Iran using real-time PCR method. A total of 71 hydatid cysts were collected from infected 
sheep, goat, and cattle slaughtered in Isfahan, Iran during 2013. DNA was extracted from protoscolices and/or germinal 
layers from each individual cyst and used as template to amplify the mitochondrial cytochrome c oxidase subunit 1 gene 
(coxl) (420 bp). Five cattle isolates out of 71 isolates were sterile and excluded from further investigation. Overall, of 66 iso- 
lates, partial sequences of the coxl gene of £ granulosus indicated the presence of genotypes G1 in 49 isolates (74.2%), 
G3 in 1 5 isolates (22.7%), and G6 in 2 isolates (3.0%) in infected intermediate hosts. Sixteen sequences of G1 genotype 
had microgenetic variants, and they were compared to the original sequence of coxl. However, isolates identified as G3 
and G6 genotypes were completely consistent with original sequences. G1 genotype in livestock was the dominant geno- 
type in Isfahan region, Iran. 
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INTRODUCTION 

Hydatid cyst, caused by Echinococcus granulosus is one of the 
most important parasitic diseases worldwide and distribute in 
many Asia countries, including Iran. This disease can cause 
high mortality in humans as well as economic losses in live- 
stock [1-3]. Nowadays, genetic studies on a partial sequence of 
mitochondrial cytochrome c oxidase subunit 1 (coxl) and nu- 
clear genetic markers suggest that E. granulosus has 10 distinct 
genotypes termed G1-G10 [4-9]. These different genotypes in- 
clude Gl and G2 as sheep strains, G3 and G5 as bovine strains, 
G4 and G6 as horse and camel strains, G7 as a pig strain, G8 
and G10 as cervid strains [10]. Recent taxonomic studies have 
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suggested E. gmnulosus in 4 distinct groups, including E. granu- 
losus sensu stricto (G1-G3 genotypes), Echinococcus ecjuinus (G4), 
Echinococcus onleppi (G5), and Echinococcus canadensis (G6-G10) 
[11]. The identification of E. granulosus strains has been carried 
out in different laboratories using various analytical methods 
(morphology, physiology, biochemistry, and molecular genet- 
ics), and all of them were proved to be useful, particularly when 
used together [12-16]. The possible susceptibility factors for 
genetic variation that occur in E. granulosus may influence many 
different phenotypic features such as the life cycle patterns, de- 
velopment rate, host specificity, geographical distribution, trans- 
mission dynamics, infectivity, antigenicity, control of diseases, 
and chemotherapeutic agents. 

PCR-based methods have been extensively used to charac- 
terize the different strains of E. granulosus [12]. On the other 
hand, various techniques such as PCR-RFLP, nested PCR and 
real-time PCR were used to determine genetic variations with- 
in E. granulosus. However, fluorescence-based quantitative PCR 
assays has not fully developed yet [17,18]. Real-time PCR is a 
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combination of PCR chemistry with fluorescent probe for de- 
tection of amplicon which shows results with curve [19]. Spe- 
cifically, EVA Green™-based real-time PCR, in association with 
melting curve analysis, can be useful for detection of nucleo- 
tide variation such as single nucleotide polymorphism (SNP). 
Different PCR products often have a different melting temper- 
ature (T ra ), and this T m is related to the percentage of guanine- 
cytosine (GC), length, and nucleotide composition [20]. In 
Iran, especially in Isfahan, sheep is the most common and im- 
portant intermediate host for E. granulosus. According to previ- 
ous studies, the prevalence of infection and cyst fertility rates 
were higher in sheep than in other slaughtered animals [21- 
23]. None of the previous studies in Iran was carried out using 
real-time PCR method for detection of E. granulosus. Therefore, 
the aim of this study was to investigate the genetic diversity of 
E. granulosus from center of Iran using real-time PCR method. 

MATERIALS AND METHODS 

Sample collection 

Seventy-one hydatid cysts were collected from infected slaugh- 
tered animals in different regions of Isfahan, including Fasa- 
ran, Khomeinishahr, and Najafabad from January 2013 until 
March 2013. Separation of hydatid cysts was performed under 
sterile condition, and then protoscolices and/or the germinal 
layer were collected from an individual hydatid cyst. Protosco- 
lices and the germinal layer were preserved in 70% (v/v) etha- 
nol and stored at -20°C until DNA extraction. 

DNA extraction 

Before DNA extraction, ethanol was removed by 3 times 
washing with sterile distilled water. Total genomic DNA was 
extracted using a genomic DNA extraction kit (Bioneer, Dae- 
jeon, Korea) according to the manufacturer's insttuctions with 
some modifications. The concentration of extracted DNA was 
evaluated by Nano Drop, and then the samples were kept at 
-20°C for further analysis. 

PCR and DNA sequencing 

Genomic DNA was amplified using specific primers for coxl 
gene, i.e., forward JB3 primer (S-TTTTTT GGG CAT CCT GAG 
GTT TAT-3) and reverse JB4.5 primer (5-TAA AGA AAG AAC 
ATA ATG AAA ATG-3'] as described by Bowles et al. [5]. PCR 
was done in a 25 pi final volume containing 10 pi master mix- 
es (Type-it HRM PCR Kit, Qiagen, Hilden, Germany), 9 pi dis- 



tilled water, 1 pi each primer, and 4 pi template DNA. 

Enzymatic reaction was performed under the following con- 
ditions: an initial denaturation step at 95°C for 10 min, follow- 
ed by 40 cycles at 95°C for 10 sec, 55°C for 30 sec, 72°C for 27 
sec, and a final extension step at 72°C for 5 min. 

The master mix provided with the kit contains the novel dou- 
ble-stranded DNA-binding fluorescent dye, Eva Green, and in- 
cludes an optimized HRM buffer, Hot StarTaq plus DNA poly- 
merase, Q-Solution, and dNTPs. Together, these components 
provided reliable PCR specificity and acceptable results, even 
with hard genomic loci. Real-time PCR was carried out in a 
Mini Opticon real-time PCR detection system (Rotor-Gene 
6,000, Hilden, Germany) under the following conditions: af- 
ter a pre-incubation step (95°C for 3 min) to activate the poly- 
merase enzyme, 40 cycles of amplification were performed, 
each one including 95°C for 10 sec, 55°C for 30 sec, and 72°C 
for 27 sec, followed by a final extension at 72°C for 5 min. The 
fluorescence signal was measured at each cycle after the exten- 
sion step. The melting experiment was performed from 55°C 
to 95°C at l°C/sec with continuous fluorescence monitoring. 
Melting peaks were visualized by plotting the first derivate aga- 
inst the melting temperature; the T m was defined as the peak 
of the curve. The quantitative detection and T m were obtained 
using the Rotor-gene 6,000 series software version 1.7 (Corbett, 
Hilden, Germany). 

Two replicates of each sample were analyzed and T m analysis 
was repeated 3 times in each run to confirm the repeatability 
of the T m assay by estimating the T m variation within PCR am- 
plification (intta-assay) and between PCR amplification (inter- 
assay), respectively. The coefficient of variation (CV) was calcu- 
lated dividing the standard deviation by the arithmetic mean 
of the measured values of T m (CV= standard deviation [SD]/ 
mean value). 

In addition, to check the uniformity of temperature of the 
cycler block, indicating the reliability of results, replicates of an 
already genotyped sample were amplified in different positions 
of the cycler block within the same amplification cycle. The in- 
tra-assay CVs represent the mean CVs for the results obtained 
from the replicates of all the E. granulosus genotypes in the same 
run. The inter-assay CVs represent the mean CVs for the results 
obtained from 3 separate runs. 

To confirm the identified genotypes, 22 samples were ran- 
domly sequenced for coxl mtDNA. Three samples that have 
been already sequenced for coxl and identified as Gl, G3, and 
G6 (accession nos. HM563013, HM563017, and HM563018, 
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respectively) were included in each PCR set as positive controls. 

Sequence homology and phylogenetic analysis 

The obtained continuous sequences were compared with 
previously published sequences of the mitochondrial coxl for 
E. granulosus genotypes in NCBI using BLAST system. Phyloge- 
netic analyses of the sequence data were inferred with maximum 
likelihood and multiple alignments and done using MEGA 5 
of the computer program CLC main workbench software (ver- 
sion 5.2.1, 2013) [24-26]. 

RESULTS 

Out of 71 isolates, 5 cattle isolates were sterile and excluded 
from further investigation. Molecular survey and PCR amplifi- 
cation were successfully performed on all 66 isolates. Of 66 
samples, 51 isolates were collected from sheep (35 livers and 
16 lungs), 8 from goat (4 livers and 4 lungs), and 7 cattle ori- 
gin (lungs). Real-time PCR using the partial sequences of coxl 
gene of 66 E. granulosus isolates showed that 49, 15, and 2 iso- 
lates were identified as Gl, G3, and G6 genotypes, respectively. 



Table 1 . Mean T m (°C) and SD calculated and intra- and inter-as- 
say CVs of the G1 , G3, and G6 genotypes of E granulosus 



£ granulosus 
genotype 


Mean 
T m CC) 


SD 


Intra-assay 
CV (%) 


Inter-assay 

CV (%) 


G1 


79.3 


0.15 


0.07 


0.12 


G3 


78.8 


0.14 


0.09 


0.14 


G6 


77.3 


0.12 


0.05 


0.12 



Sheep, goat, and cattle isolates were genotypically categorized 
as Gl (77%, 63%, and 72%), G3 (21%, 25%, and 28%), and 
G6 (2%, 12%, and 0%), respectively (data not shown). 

As shown in Table 1, real-time PCR melting curve results in- 
dicated that the mean T m of Gl, G3, and G6 genotypes were 
79.3, 78.8, and 77.3, respectively. Assessment of intra and in- 
ter-assay variability showed low and acceptable CVs (Table 1). 
Melting curve and T m of the analyzed hydatid cysts identifica- 
tion by sequencing are shown in Fig. 1. 

Sequencing results confirmed that all isolates were correctly 
differentiated by real-time PCR. All the isolates identified by 
real-time PCR were clustered along with the corresponding re- 
ference genotypes as shown in Fig. 2. The phylogenetic tree was 
divided into 2 main dades, the first dade contained 2 subclades 
corresponding to the Gl and G3 genotypes (E. granulosus sen- 
su stricto), and the second one corresponding to the isolates 
identified as G6 (£. canadensis) along with the individual refer- 
ence genotypes. Forty-nine samples with T m between 79.1°C 
and 79.5°C had identical sequences to the Glgenotype, while 
the 15 samples with T m between 78.6°C and 79.0°C had identi- 
cal sequences to the G3 genotype, and finally 2 samples with 
T m between 77.TC and 77.5°C had identical sequences to the 
G6 genotype. 

Phylogenetic analysis also showed that 16 isolates of Gl gen- 
otype had micro-variants, compared to the original sequences 
of coxl gene (accession no. M84661). Four samples identified 
as G3 and 2 samples identified as G6 were completely identical 
to the G3 and G6 reference sequences (accession nos. M84663 
and M84666, respectively). 
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Fig. 1. T m of the analyzed hydatid cysts. (A) £ granulosus G1 identified by sequencing 
and (C) £ granulosus G6 identification by sequencing. 
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Fig. 2. Molecular phylogenetic tree of 22 £ granulosus isolates of 
sheep, goat, and cow along with reference isolates based on C01 
gene sequence. The evolutionary history was inferred by using 
the Maximum Likelihood method based on the Kimura 2 param- 
eter model [19]. The tree with the highest log likelihood (-814.2241) 
is shown. Reference accession nos.: G1 , M64661 ; G3, M64663; 
G6, M84666. 

DISCUSSION 

In the present study, T m analysis was used for detection of E. 
granulosus strains (Gl, G3, and G6 genotypes) in Isfahan Prov- 
ince, Iran. Due to high sensitivity, specificity, accuracy, and re- 
peatability of this novel method, we confirmed that this tech- 
nique is useful for detection of Iranian genotypes. Because of 
0.5°C mean difference of T m between Gl and G3 genotypes, 
the Tm ranges did not overlap, therefore we could differentiate 
these 2 genotypes easily and rapidly. On the other hand, G6 
genotype could be easily distinguished from the other ones 
(Gl and G3), because of 1.5°C different T m . The accuracy of 
the T m assay was confirmed by comparing these results with 
direct sequencing. 

The advantage of using real-time PCR is that this technique 



is more cost-effective and less time-consuming than conven- 
tional PCR. Actually, with this method, a large number of sam- 
ples could be tested quickly and repeatable results could be 
achieved with significant saving in time and materials [17,27, 
28]. On the other hand, this method has some disadvantages, 
such as existing instruments with limited multiplexing capa- 
bilities, incompatibility of several systems with some fluoro- 
genic chemistries, and necessity of high technical skills for its 
launching [27]. 

Epidemiological studies showed that Gl genotype is the pre- 
dominant genotype in both human and animals [29]. Maurel- 
li et al. [17] used real-time PCR for identification of Gl and 
G2/G3 genotypes of E. gmnulosus in Italy. They found that T m 
analysis can be used for discriminating between Gl and G2/ 
G3 genotypes of E. granulosus sensu stricto. Our results were 
compatible with their findings. 

Recently, Boubaker et al. [30] in Switzerland developed a 
single-tube multiplex PCR (mPCR) for discrimination of the 3 
levels of the E. g-anulosus complex genotyping. They found E. 
gi-anulosus sensu stricto (G1/G2/G3), E. granulosus eauxnus (G4), 
and E. canadensis (G7) as predominant genotypes in their study. 
This study showed that this method could provide the unique 
opportunity to address directly speciation and genotyping with- 
in the framework of large-scale studies. 

Recent studies conducted analysis of various strains in the 
Middle East. For example, a study by Adwan et al. [12] in Pal- 
estine showed the presence of G1-G3 genotypes in this region 
using coxl loci. In another study, Simsek et al. [31] reported a 
G1/G3 cluster (£. granulosus sensu stricto) in Turkey. 

In recent years, several studies reported different genotypes 
of E. granulosus from intermediate and definitive hosts, includ- 
ing the sheep, goat, cattle, camel, buffalo, human, and dog in 
different parts of Iran [4,22,32,33]. The results of this study 
were similar to previous studies in Iran showing that genotypes 
Gl, G3, and G6 are present in intermediate hosts [29]. In the 
few studies, genotypes Gl, G3, and G6 were simultaneously 
reported [4,34]. Recently, Rostami et al. [29] used the post real- 
time PCR with HRM method and differentiated Gl, G3, and 
G6 genotypes of sensu lato in Iran. Rajabloo et al. [35] report- 
ed Gl and G6 genotypes from goats in Iran. A similar study in 
Ham Province on hydatid cysts isolated from animals and hu- 
mans showed E. granulosus sensu stricto (G1-G3) in this region 
[36], A study on 120 hydatid cysts isolated from goats in Ma- 
zandaran Province, Northern Iran, by Youssefi et al. [37] show- 
ed the Gl genotype as the dominant strain. Among various 
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studies, Rostami et al. [29] reported Gl, G3, and G6 genotypes 
in the intermediate hosts that are similar to our results. 

In conclusion, we believe that real-time PCR that has been 
developed in the present study, provides powerful means for 
molecular and epidemiological studies on human and animal 
E. granulosus infections. Therefore, we suggest that the extensive 
use of this technique with high sample size can help to find 
undetected genotypes using different molecular targets. 
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